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I .  INTRODUCTION 


The  loading  on  equipment  and  structures  caused  by  shock  diffraction 
and  drag  have  been  of  particular  interest  to  the  Army  ever  since  the  in¬ 
troduction  of  nuclear  weapons.  Because  of  the  inherent  three-dimensional 
(3-D)  nature  of  such  loading,  the  more  commonly  accepted  method  for 
ascertaining  this  loading  has  been  through  field  testing,  dating  from 
atmospheric  nuclear  tests  to  current  tests  of  large-scale  high  explosives. 
Because  of  the  high  cost  and  necessarily  complex  planning  involved,  these 
tests  are  held  infrequently.  Unpredictable  anomalies  such  as  jetting 
frequently  occur,  often  with  the  intended  targets  receiving  loadings 
considerably  different  from  that  which  had  been  expected.  Space  for 
targets  and  data  acquisition  equipment  is  invariably  at  a  premium 
during  these  tests,  leaving  potentially  valuable  tests  untried,  and 
worthwhile  data  ungathered. 

As  pointed  out  by  Taylor1  in  1972,  the  accepted  methods  for 
computing  shock  diffraction  loading  on  simple  structures  as  outlined  in 
the  US  Army  Technical  Manual  TM-S-856-12,  hereinafter  referred  to  as 
the  Technical  Manual,  are  no  longer  sufficiently  accurate  for  studies 
in  blast  hardening  and  vulnerability  assessment.  The  models  outlined 
in  the  Technical  Manual  treat  a  given  3-D  structure,  such  as  a  rectan¬ 
gular  parallelepiped,  as  being  locally  two-dimensional.  Loading  is  then 
computed  based  on  the  speed  at  which  a  rarefaction  wave  travels  across  a 
characteristic  dimension  of  the  given  target  face.  During  the  last 
several  years,  numerical  methods  in  the  form  of  hydrodynamic  computei 
codes  (hydrocodcs)  for  simulating  3-D  shock  diffraction  over  obstacles 
have  boon  advanced  to  the  point  where  they  can  be  used  to  complement 
and,  in  some  cases,  replace  the  use  of  experiments  and  models  such  as 
those  in  the  Technical  Manual. 

The  objective  of  the  present  work  is  to  illustrate  the  capabilities 
of  two  hydrocodes  for  predicting  shock  diffraction  loading  and  to  com¬ 
pare  them  with  two  semiempirical  models.  The  hydrocodes  used  for  the 
comparisons  arc  HUM.  and  BAAL.  One  semiempirical  model  is  that  from 
the  Technical  Manual,  and  the  other  is  an  improved  model  suggested  by 
Taylor1 . 

A  reference  problem  is  defined:  the  S-280  Llectrical  Lquipment 
Shelter  struck  by  a  34. S  kPa  (5.0  psi)  overpressure  step  shock.  Im¬ 
provements  to  the  models  for  determining  loading  in  the  Technical 


Part 


J{/.  ,/.  Taylor,  "A  Method  for  Mediating  Blast  Iroada  !hanng  the 
H? fraction  Pham,"  The  Shock  and  Vibration  Bulletin,  HP- 
./  of  6,  Shook  and  Viln'ation  Center,  tiaoal  Research  Lalwmtory, 
Washington,  DC,  p.  135  (Jan  ?3). 

S  "Design  of  Structures  to  Rosie t  the  Rf foots  of  /!< ante  Wcaj p* 
Amy  Corps  of  (Engineers  EN-lllO-m-41*,  !  duly  2959,  (VllCLA^msD) . 
(Republished  as  TMS-SSC-l  in  19B5.) 


HtscjuuiKt  aunc-iior  filmed 


Manual  are  suggested;  these  improved  models  give  good  agreement  with 
the  hydrocou.  results.  The  BAAL  hydrocode  has  previously  been  shown 
to  compare  well  with  experimental  shock  loading  on  a  similar  rectan¬ 
gular  parallelepiped  at  the  same  shock  strength.  Finally,  a  discussion 
is  presented  on  the  net  rotational  moments  on  the  shelter  as  predicted 
by  the  several  methods. 


II.  REFERENCE  PROBLEM 


A.  Target 


The  target  is  an  S-280  Electrical  Equipment  Shelter  (hereinafter 
referred  to  as  the  shelter),  3.62  m  wide,  2.17  m  deep,  and  2.11  m  high. 
The  shelter  is  sitting  on  the  ground  with  one  of  the  larger  of  the  two 
different  sized  faces,  defined  here  as  the  front  face,  oriented  so  that 
it  is  normal  to  the  velocity  vector  of  an  oncoming  one-dimensional  (l-P) 
shock  wave. 

B.  Shock  and  Ambient  Conditions 


The  shock  wave  is  a  step  shock  with  an  overpressure  p*  =  34.5  kPa 
(5.0  psi).  Ambient  conditions  ahead  of  the  shock  are  temperature 
T1  =  288. 2°K  (518. 7°R),  pressure  Pj  =  101.3  kPa  (14.7  psi),  and 

density  Pj  =  1.225  kg/ra3  (.0765  lbra/ft5).  The  air  both  ahead  of  and 

behind  the  shock  is  assumed  to  bo  a  polytropie  gas,  with  a  ratio  of 
specific  heats  y  =  1.4.  Throughout  the  remainder  of  this  report,  the 
interaction  of  this  shock  wave  with  the  shelter  described  above  will 
frequently  ho  referred  to  as  ''the  present  case." 

C.  One-Dimonsional  Chock 

A  simplified,  locally  1-0  model  that  can  be  used  to  represent  the 
flow  field  at  the  time  that  the  incident  shock  reflects  from  the 
shelter's  front  face  is  shown  u»  Figure  1.  The  flow  field  shown  is 
for  a  shock  tube  with  a  closed  end.  Region  1  represents  the  ambient 
air;  region  2,  the  3ir  behind  the  incident  shock  wave;  region  3,  the 
shock  tube  reservoir;  and  region  4,  the  accelerated  flow  region  between 
the  expansion  wave  and  the  contact  surface.  Finally,  region  5  is  the 
region  behind  the  shock  wave  which  reflected  off  the  rigid  wall  that 
represents  the  shelter's  front  face.  The  subscripts  used  with  flow 
field  variables  in  this  report  indicate  the  particular  region  of 
interest . 

The  values  of  the  flow  field  variables  for  the  1-t)  check  can  be 
calculated  from  the  shock  tube  relations  with  the  shock  and  ambient 
conditions  specif icd  above.  These  1-1)  variables  arc  summari zed  in  the 
table  in  Appendix  A.  The  theoretical  value  of  the  peak  reflected 
overpressure  of  78.5  kPa  calculated  with  the  1-t)  model  is  an  excellent 
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Figure  l.  Sieplificd  onc-dimensional  flow  field. 
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reference  point  to  use  in  assessing  the  validity  of  the  hydrocode 
solutions  of  the  real  3-D  case.  The  initial  reflection  process  at  any 
given  point  on  the  front  face  is  1-D  until  the  expansion  waves  coming 
in  from  the  top  and  side  edges  of  the  front  face  arrive  at  that  point. 


A.  Background 


III.  SKMI EMPIRICAL  MODF.I.S 


Prior  to  presenting  the  results  of  the  hydrocode  computations,  it 
is  worthwhile  to  discuss  the  existing  semiempirical  models  used  to. 
determine  shock  loading  on  3-D  structures.  This  allows  not  only  a 
comparison  between  the  loading  predicted  by  the  models,  but  also  allows 
later  comparison  to  the  loading  predicted  by  the  hydrocodes.  The  chosen 
semiempirical  models  are  (1)  those  from  the  Technical  Manual’,  and  (2) 
the  improved  models  suggested  by  Taylor1. 

B.  Technical  Manual 


The  equations  presented  here  are  for  a  step  shock,  and  any  mention 
of  the  terms  "shock11  or  "shock  wave"  in  the  following  discussion  arc  to 
be  construed  as  referring  to  a  step  shock  unless  specifically  noted  to 
be  otherwise.  The  semiempirical  models  in  the  Technical  Manual  are  for 
decaying  shock  waves,  but  the  relations  can  be  easily  modified  to  relate 
to  step  shocks.  The  derivations  of  the  step  shock  equations  from  the 
decaying  shock  wave  equations  in  the  Technical  Manual  are  given  in 
Appendix  ft. 

The  Technical  Manual  models  of  interest  here  assume  that  the  shock 
wave  strikes  a  rectangular  parallelepiped  with  the  velocity  vector  of 
the  shock  wave  normal  to  one  of  the  faces  of  the  parallelepiped.  That 
face  is  defined  as  the  front  face.  A  semiempirical  model  which  shows 
the  shock  reflection  and  subsequent  expansion  wave  interact i or, s  is 
given  for  the  front  face.  A  semi empirical  model  is  given  for  the  back 
face  which  shows  the  shock  expansion  around  the  po-degree  comers  with 
the  top  and  side  faces.  A  combined  model  is  given  for  the  top  and  side 
faces,  showing  shock  movement  along  those  faces  plus  wave  interactions 
from  the  front  and  back  faces. 

These  models  from  the  Technical  Manual,  modified  for  step  shocks, 
are  discussed  below.  The  predictions  made  using  these  models  are  also 
discussed  and  presented  for  later  comparison  with  other  models  and  with 
the  h/drocodc  results.  Any  numerical  values  such  as  sound  speed,  shock 
speed,  or  overpressure  refer  specifically  to  the  interaction  of  the 
reference  34.5  kPa  Step  shock  with  the  snelter  (the  present  case). 

1.  Front  Face.  The  model  in  the  Technical  Manual  for  calculating 
the  load ing  un~The~~ front  face  is  based  on  an  assumed  clearing  time  tf, 
where 


CD 


t 

c  c5 

Here,  is  the  sound  speed  on  the  front  face  of  the  shelter  after  the 

shock  has  reflected  from  it;  for  the  present  case,  c^  =  369.6  m/s.  The 

value  to  be  used  for  h  is  the  smaller  of  either  the  height  of  the  shelter 
or  one-half  the  width;  for  the  present  case,  h  =  1.81  m.  This  gives  a 
clearing  time  for  the  shelter  t  =  14.69  ms,  where  t  is  "...  the  time 

required  to  clear  the  front  wall  of  reflection  effects."2  During  this 

time,  the  average  overpressure  on  the  front  face,  p^,  decreases  from  the 

peak  reflected  overpressure  to  a  value 

p*  =  (p2  -  px)  +  0.85  ~  P2u22‘  ^ 

Here,  p2  is  the  pressure  behind  the  shock,  p2  is  the  density  behind  the 

shock,  and  u2  is  the  particle  velocity  behind  the  shock.  For  the 

present  case,  the  peak  reflected  overpressure  on  the  front  face  is 

(p5  -  Pl)  or  78.5  kPa  (11.4  psi) ,  and  p*  is  37.9  kPa  (5.49  psi);  p$  is 

the  pressure  behind  the  reflected  shock.  Figure  2  shows  the  history  of 
the  average  overpressure  on  the  front  face  of  the  shelter  calculated 
using  the  model  in  the  Technical  Manual.  The  term  "average  overpressure" 
refers  to  a  spatial  average  on  the  shelter  face  at  a  given  point  in 
time,  or  at  a  series  of  points  in  time.  In  the  later  discussions 
concerning  the  hydrocode  computations,  average  overpressure  means  a 
weighted  spatial  average,  with  the  weighting  based  on  the  areas  of  the 
flow-field  cells  adjacent  to  the  shelter  surface.  The  other  data 
shown  in  Figure  2  for  the  front  face  will  be  discussed  later. 

2.  Back  Face.  The  time  required  for  the  incident  shock  wave  to 
arrive  at  the  plane  of  the  back  wall  is 

t  =  ~  (D 

cd  W.  ’ 

i 

where  D  is  the  shelter  depth,  and  W.  is  the  velocity  of  the  incident 
shock.  For  the  present  case,  D  =  2.17  m  and  =  386.8  m/s,  so 
t.  =  5.616  ms.  The  rise  time2  required  for  the  average  pressure  on  the 
back  face  to  go  from  ambient  at  time  t^  to  its  peak  value  is 


(t  j  =  — 
^  bJnse  Cj 
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where  a  =  4.0,  is  the  sound  speed  ahead  of  the  shock,  and  h  is 

defined  as  before.  Here,  h  =  1.81  m  and  c,  =  340.3  m/s,  so  (t.)  .  - 

1  '  b  rise 

21.28  ms.  The  peak  average  overpressure7  on  the  back  face  is 


where 


^b^max 


Pi) 


1 


+ 


8 


(s) 


(0) 


or  (p£)  =  29.3  kPa. 

vrVmax 

Because  the  shock  wave  considered  here  is  a  step  shock,  the  model 
in  the  Technical  Manual  implies  that  the  average  overpressure  on  the 
back  fane  would  remain  at  (?pmax  indefinitely.  The  overpressure 

history  calculated  using  this  model  is  represented  by  the  solid  line 
labeled  "BACK"  in  Figure  2. 

Figure  3  shows  the  history  of  the  average  pressure  difference 
between  the  front  and  back  faces  of  the  S-280  shelter  using  the 
Technical  Manual  model.  Tho  other  data  in  Figures  2  and  3  will  be 
discussed  later. 

3.  Top  and  Side  Faces.  From  the  time  t  =  0.0  s,  when  the  incident 
shock  wave  r ea ches  the  plane  of  the  front  face  of  the  target,  to  the 
time  t^  at  which  t.ho  shock  wave  reaches  the  plane  of  the  back  face,  the 

average  overpressure  on  tho  top  face  varies  linearly7  from  aero  to 


=  (ps  •  pi) 


P-*  *  Pi 

0.9  ♦  0.1  |  1.0  -  ----- 

Pj 


(7) 


For  the  present  case,  p*  *  32.5  kl’n  at.  t(j  »  5.616  ms.  To  account  for 

vortex  growth,  shedding,  and  subsequent  travel  down  the  top  face,  the 
model  predicts  a  local  minimum  in  the  average  overpressure  at  the  time 
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which  for  the  present  case  is  28.08  ms.  At  that  time,  the  average 
overpressure  at  the  top  face  is  calculated2  as  either 


or 


0.5  +  0.125  2 


P;  -  PiV 

pi  ) 


(9.2) 


whichever  is  the  lesser  of  the  two.  At  this  time,  p*  =  25.2  kPa  for 

1  top 

the  present  case.  After  this  time,  there  is  a  linear  rise  in  the 
average  overpressure  to  a  value 


from  the  time  t 

pmxn 


For  the  present  case 
shock  overpressure. 


Ptop  =  ^2  ‘  h 
to  a  time 


t  -  t  iiD 

pmin  +  W. 

t  =  98.28  ms  and  p* 

1  top 


(10) 


(ill 

34.5  kPa,  the  incident 


According  to  the  method  outlined  in  the  Technical  Manual,  the 
geometry  of  the  shelter  is  such  that  the  times  and  average  overpressures 
for  the  top  face  also  apply  directly  to  the  side  faces  of  the  shelter. 
Figures  4  and  5  show  the  overpressure  history  for  the  top  and  side 
faecs,  respectively,  of  the  shelter  for  the  present  case  calculated 
with  the  above  method.  The  hydrocode  data  and  the  proposed  model  shown 
in  those  figures  will  be  discussed  later. 


C.  Taylor's  Modct 

1.  Front  Face.  Taylor*  suggested  an  alternate  model  to  that  of 
the  TechnI caTlfanua 1  to  estimate  the  loading  on  the  front  face  of  an 
obstacle,  based  on  the  number  of  rarefaction  wave  crossings.  For  a 
3-P  object  such  as  the  shelter,  there  will  be  a  succession  of  rarefac¬ 
tion  waves  originating  from  the  top  of  the  front  face  and  reflecting 
from  the  bottom  boundary,  plus  a  succession  of  rarefaction  waves 
originating  from  the  side  of  the  front  face  and  reflecting  from  the 
symmetry  plane.  The  reflected  waves  will  also  be  rarefaction  waves. 

In  an  actual  occurrence  of  loading  through  shock  diffraction,  there 
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Figure  4.  Top  face  average  overpressure. 
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Figure  5.  Side  face  average  overpressure. 
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would  be  considerable  interaction  between  the  crossing  expansion  waves 
Because  Taylor  intended  that  the  model  be  a  simple  approximation,  the 
waves  are  assumed  not  to  interact  with  one  another.  As  a  further  sim¬ 
plification,  the  rarefaction  wave  speeds  are  assumed  to  be  equal  to  the 
sound  speed  in  the  reflection  region  immediately  after  the  incident 
shock  reflects  from  the  front  face. 

Using  Taylor's  model,  the  required  crossing  time  for  an  expansion 
wave  running  from  the  top  to  the  bottom  of  the  front  face  is  5.704  ms; 
for  an  expansion  wave  running  between  the  side  edge  and  the  symmetry 
plane  the  crossing  time  is  4,897  ms.  These  times  are  based  on  a  shelter 
height  of  2.11  m,  a  one-half  width  of  1.81  m,  and  a  sound  speed 
c5  =  m/s.  For  each  expansion  wave  crossing,  the  average  over¬ 

pressure  on  the  front  face  is  then  calculated  as  a  percentage  of  the 
initial  reflected  overpressure  on  the  front  face.  The  values  for  average 
overpressure  using  Taylor's  fitted  curve1  are  summarized  in  Table  1  and 
shown  in  Figure  2. 


TAl ’ 7  1.  FRONT  FACE  AVERAGE  OVERPRESSURE,  TAYLOR'S  MODEL 


Time 

(ms) 

Rarefaction 
Wave  Number* 

Wave 

Description 

%  Reflected 
Overpressure 

Overpressure 

(kPa) 

0.0 

— 

- 

100.0 

78.5 

4.9 

1 

Side  to  Symmetry 
Plano 

86.9 

68.3 

5.7 

2 

Top  to  Bottom 

68.5 

55.8 

9.8 

3 

Symmetry  Plano 
to  Side 

57,3 

45.0 

11.4 

4 

Bottom  to  Top 

52.8 

41.5 

14.7 

5 

Side  to  Symmetry 
Plano 

SO. 4 

39.6 

^ Model  stops  at  crossing  number  5. 

Taylor's  model  indicates  a  more  rapid  unloading  of  the  front  face 
during  the  diffraction  phase  than  does  the  Technical  Manual.  (The 
diffraction  phase  on  the  front  face  lasts  from  the  initial  shock  arrival 
until  all  wave  interactions  on  that  face  have  ended.)  The  two  models 
are  nearly  the  same  for  the  drag  phase  on  the  front  face  (t  >  15  ms). 
(The  drag  phase  begins  nt  the  end  of  the  diffraction  phase  for  a  given 
face.)  The  diffraction  phase  for  the  whole  shelter  ends  when  there 
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are  no  longer  any  wave  interactions  anywhere  on  the  shelter;  this  ends 
at  t  =  20  ms.  This  and  the  hydrocode  data  shown  in  Figure  2  will 
be  discussed  later. 

2.  Back  Face.  Taylor's1  data  indicate  that  the  constant  a  used 
in  Equation  (4)  should  be  changed  to  2.5  from  4.0.  This  gives  a  more 
rapid  pressure  rise  time  of  13.30  ms  on  the  back  face  as  compared  with 
that  calculated  with  the  model  in  the  Technical  Manual .  The  overpressure 
history  for  the  back  face  of  the  shelter  calculated  with  Taylor's 
method  is  also  shown  in  Figure  2.  It  shows  much  more  rapid  loading  of 
the  back  face  than  does  the  Technical  Manual  model;  the  late-time  loading 
is  the  same.  Figure  3  shows  the  history  of  the  average  pressure 
difference  between  the  front  and  back  faces  of  the  shelter  using  the 
models  outlined  by  Taylor.  The  net  loading  predicted  by  Taylor's  models 
is  much  less  than  that  for  the  Technical  Manual  for  nearly  all  of  the 
diffraction  phase.  The  hydrocode  data  shown  in  Figure  3  will  be 
discussed  in  a  later  section. 

3.  Top  and  Side  Faces.  Taylor  did  not  suggest  an  improved  method 
for  computing  the  average  overpressure  on  the  top  and  side  faces  of  an 
object  such  as  the  shelter. 


IV.  HYDRQCQDES 


A,  Background 

Until  relatively  recently,  models  such  as  those  outlined  in  the 
Technical  Manual  were  the  accepted  means  for  calculating  the  loading 
on  a  structure  being  struck  by  a  shock  wave.  Taylor1  pointed  out  the 
more  serious  deficiencies  in  the  simple  models  being  used,  particularly 
the  poor  modeling  of  three-dimensional  effects,  and  suggested  a  simple 
but  effective  model  based  on  wave  interactions.  A  more  generalised 
version  of  Taylor's  wave  interaction  model  has  since  been  adapted  in 
the  TKUCK  code  by  Hobbs,  et  al.3  Two  features  of  this  model  are  the 
ability  to  model  oblique  loading  and  more  complex  target  shapes. 

However,  it  was  apparent  that  a  more  general  and  accurate  computa¬ 
tional  model  titan  that  offered  by  the  simple  models  was  needed.  At  *hc 
time  of  Taylor's1  paper,  hydrodynamic  computer  codes  capable  of  solving 
the  Kuler  equations  (and  in  limited  cases  the  Navier -Stokes  equations) 
for  two-dimensional  flow  were  fairly  well  established.  Advances  were 
also  being  made  in  the  development  of  hydrodynamic  computer  codes 


3tt.  P.  Hobbs,  J.  P.  Walsh,  G.  Tartarian,  W.  N.  Lev,  and  V.  Wu,  M?&rK  - 
A  Digital  Cemjmter  Program  foe  Calculating  the  Respmee  of  Ar**g 
Vehicles  to  Blast  Waves,  *  Ballistic  Research  Laboratertt  ConUxtct 
Report  Do.  00SG3,  ApHi  W?8,  tUBCLASSIPllCD) .  (At)  «AOSn5J7) 
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fXrafnw  8  three-dimensional  problems.  The  Ballistic  Research 

Laboratory  (BRL I  wanted  to  obtain  a  reliable  3-D  hydrocode  for  shock 
diffraction  loading  and  free  air  blast  calculations.  The  two  hydrocodes 
tested  m  detail  are  described  below. 


B.  BAAL 


One  of  the  several  hydrodynamic  computer  codes  being  developed  at 
the  Los  Alamos  Scientific  Laboratory  (LASL)  was  the  implicit,  arbitrary- 
Lagrangian-Eulerian  (ALE)  hydrocode  BAAL4.  The  computing  methods  used 
m  BAAL  are  described  for  two-dimensional  flow  by  Hirt,  Amsden,  and 
Cook  ,  and  also  by  Amsden  and  Hirt6.  The  extension  of  these  methods  to 
three-dimensional  flow  is  described  by  Pracht7.  As  a  result  of  dis¬ 
cussions  with  representatives  of  the  Ballistic  Research  Laboratory 
(BRL),  LASL  was  contracted  to  perform  a  computation  using  the  BAAL 
hydrocode,  simulating  the  three-dimensional  shock  diffraction  loading 
experiment  performed  by  Taylor1  at  the  BRL  using  a  solid  rectangular 
parallelepiped  as  a  target.  The  results  of  this  computation  have  been 
reported0 .  These  results  have  also  been  reported  in  the  open  literature 
along  with  a  discussion  of  the  solution  technique9.  The  BAAL  calculation 
involved  a  solution  of  the  inviscid  F.uler  equations,  with  artificial 
viscosity  used  to  stabilize  flow  in  regions  of  deceleration.  (The  BAAL 


S.  I*raaht  and  J.  V,  Braakbill ,  "BAAL:  A  Cexle  for  Calculating 
Three-Vimn&ional  Fluid  Flew  at  All  Speeds  m'th  an  Fulerian- 
tog  rang  ion  Computing  Meek, "  toe  Almas  Soimtifia  Laboratory 
Report  IA-6S42,  August  1876. 

C.  M.  Hirt,  A,  A.  Amsden,  and  J.  L.  Cook,  "An  Arbitrary  Ligmngian- 
Sulerian  Computing  Method  for  All  Floo  Speeds , "  A.  Cmp.  Vbnt.  H, 
1374,  pp.  9S7-S63. 

A.  A.  Amsden  and  C,  W.  Hirt,  "TAQVI:  An  Arbitrary  tograngian- 
Sulerian  Computer  Pi'ogrm  for  Fluid  Ftm»  at  All  Speeds,"  Loo  A  tome 
Sai-mtifia  Laboratory  Report  lA-bWO,  March  13??. 


w.  H,  ftnoht,  "Calculating  Tkree-tHmntiioml  Fluid  Flew  at 
oith  an  Fu lev utn- Laffremgian  Cittreuting  Mesh,"  A.  Comp  Phya. 
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code  has  the  capability  of  solving  the  full  Navier-Stokes  equations  for 
some  flow  situations.)  In  general,  the  BAAL  computation  showed 
excellent  agreement  with  the  experimental  pressure  data  for  both  the 
front  and  back  faces.  Experimental  data  were  not  taken  on  the  top  and 
side  faces.  As  a  result  of  the  confidence  gained  in  the  BAAL  code 
from  this  comparison,  LASL  was  contracted  to  perform  a  second  computa¬ 
tion.  This  was  for  an  S-280  Electrical  Equipment  Shelter  struck  by  a 
34.5  kPa  (5.0  psi)  step  shock  wave,  the  reference  case  discussed 
earlier.  The  results  of  this  computation  have  been  reported  previously10. 
Using  the  data  generated  in  the  BAAL  computation10  for  the  34.5  kPa 
overpressure  shock  wave  diffracting  over  the  shelter,  and  using  a  simple 
scaling  model  suggested  by  the  BRL  to  extrapolate  the  34.5  kPa  data  to 
estimate  pressure  histories  for  other  comparable  shock  overpressures, 
Calligeros  et  al . 1 1  predicted  responses  of  the  shelter  to  various  shock 
waves  for  Event  Dice  Throw.  Ethridge12  has  used  Taylor's  results1 
together  with  other  experimental  data  to  construct  correlation  functions 
for  blast  diffraction  loading  on  the  front  and  rear  surfaces  of  a  rec¬ 
tangular  parallelepiped.  Those  correlation  functions  are  compared12 
with  the  data  of  Taylor1  and  Lottero10,  showing  good  agreement. 

C.  HULL 


During  approximately  the  same  period  that  the  three-dimensional 
BAAL  computations  were  being  run,  a  three-dimensional  capability  was 
being  added  to  the  HULL1 3  hydrodynamic  computer  code  at  the  US  Air 
Force  Weapons  Laboratory  (AFWL) .  The  major  differences  between  the  HULL 
and  BAAL  hydrocodes  are  that  HULL  solves  only  the  inviscid  Ruler 
equations,  uses  an  explicit  finite  difference  algorithm,  and  is  re¬ 
stricted  to  a  fully  Eulerian  computing  mesh  comprised  of  cells  which 
are  rectangular  parallelepipeds.  The  AFWL  wished  to  verify  the  new 


R.  E.  Lottero ,  " Computational  Predictions  of  Shook  Diffraction  Loading 
on  an  S-280  Electrical  Equipment  Shelter ,"  Ballistic  Research 
laboratory  Memorandum  Report  No.  2599,  March  1976 ,  (UNCLASSIFIED) . 

J.  P.  Walsh ,  and  R.  P.  Yeghiayan ,  "Structural 
Modeling  and  Response  of  Comand ,  Control  and  Communication  Shelter 
Syotema  for  Event  Dice  Throw, "  Ballistic  Research  lat'oratorij 
Contract  Report  No.  00421 ,  March  1980 ,  (UNCLASSIFIED) .  ///V085759) 

J21V.  R.  Ethridge ,  "Blast  Diffraction  Loading  on  the  Front  and  Rear 
Surfaces  of  a  Rectangular  Parallelepiped,"  Ballistic  Research 
Laboratory  Memorandum  Report  No.  2784,  September  1977,  (UNCLASSIFIED) . 

(AD  #A05U312) 

liM.  A.  Fry,  RSE.  Durret,  G.  P.  Ganong,  D.  A.  Matuska ,  M.  P.  S tucker , 

B.  S.  Chambers,  C.  E.  Needham,  and  C.  D.  Westmoreland ,  "The  NULL 
Hydrodynamics  Computer  Code,"  US  Air  Force  Weapons  laboratory 
Technical  Report  76-183,  September  19??,  (UNCLASSIFIED) . 
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three-dimensional  HULL  against  existing  experimental  data  or  proven 
codes.  The  BRL  wished  to  gain  access  to  another  reliable  three-dimen¬ 
sional  hydrocode,  because  it  appeared  that  BAAL  might  not  be  immediately 
available  for  production  work.  The  AFWL  supplied  the  computing 
facilities  and  the  three-dimensional  HULL  code;  the  BRL  supplied  the 
test  problem  and  the  analysis  of  the  hydrocode  output. 

The  problem  chosen  for  the  HULL  computation  was  the  reference 
problem  described  earlier,  identical  to  that  run  in  the  second  BAAL 
computation10  including  the  flow  field  grid  formulation.  This  was 
chosen  to  gain  corroboration  for  the  BAAL  computation  for  the  shelter 
in  time  for  Event  Dice  Throw. 


V.  HYDROCODE  COMPUTATIONS 

Before  discussing  the  results  of  the  hydrocode  calculations,  the 
general  nature  of  the  computations  will  be  described.  This  section  and 
its  associated  Appendix  C  together  describe  the  direction  and  index 
conventions,  the  computational  grid,  the  computation  identifications, 
various  flow  parameters,  computer  mainframes  and  central  processor 
unit  (CPU)  times,  initial  shock  locations,  and  other  such  characterizing 
data  concerning  the  computations. 

Because  the  HULL  and  BAAL  codes  use  different  conventions  for 
naming  directions  and  indices,  the  BAAL  results  and  the  flow  field 
diagram  are  reproduced  in  part  using  the  HULL  hydrocode  conventions. 
Figure  6  shows  a  3-D  view  of  the  shelter  in  the  computational  flow 
field.  The  computational  grid  is  discussed  in  detail  in  Appendix  C. 

While  there  is  just  one  BAAL  computation  for  the  shelter,  there 
are  two  complete  HULL  computations,  designated  HULL  19.8063  (hereafter 
referred  to  as  HULL  A)  and  HULL  19.8067  (hereafter  referred  to  as 
HULL  B) .  There  is  also  an  earlier  incomplete  HULL  computation, 
designated  HULL  202.1  (hereafter  referred  to  as  HULL  C) ,  that  contains 
some  interesting  features.  All  three  HULL  runs  are  for  the  same 
34.5  kPa  step  shock  wave.  The  initial  formulations  of  the  four 
computations  are  summarized  in  Table  2. 

Analysis  by  the  BRL  revealed  that  the  computation  for  HULL  C  was 
not  valid  after  19  ms  simulated  problem  time  because  of  the  arrival  of 
an  artificial  wave  at  the  front  face.  The  reflected  wave  from  the 
front  face  had  traveled  upstream,  eventually  arriving  at  the  upstream 
transmissive  boundary  through  which  the  HULL  code  had  been  feeding  the 
stop  shock.  Because  the  algorithm  for  the  upstream  transmissive 
boundary  was  formulated  to  function  as  a  simple  positive  image 
boundary,  the  initially  constant  input  wave  was  modified  by  the  arriving 
compression  wave  and  its  following  expansion  wave.  Ultimately,  this 
artificial  interaction  at  the  upstream  transmission  boundary  caused  the 
incoming,  originally  steady,  wave  to  be  almost  completely  turned  off, 
with  the  artificial  wave  arriving  at  the  front  face  by  19  ms  simulated 
problem  time. 


f-'tgarr-  6,  S-2S0  electrical  equipment  Shelter  in  the 
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With  the  information  from  the  anal-  his  of  HULL  C  and  from 
extensive  analyses  of  other  problems  b-  AFWL,  improvements  to  3-D  HULL 
were  made  by  AFWL.  The  problems  labeled  HULL  A  and  HULL  B  were  run  at 
AFWL  in  a  cooperative  effort  by  the  BRL  and  AFWL  using  the  improved 
3-D  HULL.  These  two  HULL  runs  are  identical  except  that  for  HULL  A  the 
shock  is  started  at  seven  flow  field  cells  upstream  of  the  S-280 
shelter  front  face,  whereas  for  HULL  B  the  shock  is  started  at  the 
front  face  of  the  shelter,  as  was  the  case  for  the  BAAL  run  and  HULL  C. 
Artificial  viscosity  was  used  for  both  the  BAAL  run  and  HULL  C;  no 
artificial  viscosity  was  used  in  HULL  A  or  HULL  B.  Also,  HULL  C  was 
run  with  a  Courant-Friedrichs-Lewy  (CFL)  stability  factor  n  =  0.15, 
while  both  HULL  A  and  HULL  B  used  r,  =  0.50. 


VI.  LOADING  ON  SHELTER  FACES 

A.  Remarks 

This  section  contains  detailed  discussions  of  the  loading  on  the 
shelter,  with  individual  sections  on  the  front,  back,  top,  and  side 
faces.  The  hydrocode  results  and  the  predictions  made  from  the 
Technical  Manual's  ~.-d  Taylor's  semiempirical  models  are  compared  with 
one  another. 

B.  Front  Face 

Figure  7  shows  a  comparison  of  the  average  overpressure  on  the 
front  face  of  the  shelter  computed  with  BAAL;  HULL  runs  A,  B,  and  C; 
the  Technical  Manual;  and  Taylor's  model.  The  two  simple  models  both 
show  the  theoretical  peak  overpressure  of  78.5  kPa  at.  time  zero.  The 
first  data  point  for  the  BAAL  computation  is  at  t  =  1.S7  ms,  showing 
an  average  overpressure  of  74,3  kPa  for  the  front  face,  5,4  percent 
below  the  theoretical  peak  value.  No  pressure  data  prior  to  this  time 
were  furnished  to  the  BRL,  and  hence  the  character  cf  the  BAAL  computa¬ 
tion  from  0  <  t  <1.571  ms  is  unknown. 

HULL  C  is  the  HULL  computation  that  most  closely  resembles  the 
BAAL  computation  in  its  general  nature.  These  two  computations  use 
essentially  the  same  grid  (See  Appendix  C);  they  are  both  started  with 
the  shock  at  the  front  face.  Both  computations  use  a  form  of  artificial 
viscosity,  applying  it  only  in  rcgicrs  of  deceleration,  although  the 
actual  form  of  the  BAAL  viscosity  function4'  is  different  from  that  used 
in  the  HULL  code14.  As  can  be  seen  i;.  Figure  7,  the  HULL  C  computation 
appears  to  be  stable  from  time  zero  up  to  19  ms,  after  which  the 

14 P,  w.  hum ,  H.  J.  Happ  III ,  and  C.  E.  Needham ,  "Development  of  an 
Artificial  Viscosity  Function US  Air  Force  Weapon ft  Moratory 
Technical  Beport  77-53,  September  19?? ,  (UNCLASSIFIED) . 
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solution  drops  off  significantly  due  to  the  arrival  of  the  artificial 
wave  discussed  earlier.  The  computation  after  19  ms  is  not  plotted. 

The  HULL  C  computation  reaches  a  peak  average  overpressure  of  approxi¬ 
mately  69.0  kPa  at  1.3  ms,  12.1  percent  less  than  the  theoretical  peak. 
From  a  time  of  2.5  ms  to  17  ms,  the  BAAL  and  HULL  C  computations  are 
nearly  identical.  After  this  time,  the  BAAL  computation  tends  to  rise 
to  values  above  the  others,  whereas  the  HULL  C  computation  begins  to 
converge  with  the  other  HULL  computations  until  it  is  destroyed  by  the 
artificial  wave. 

The  agreement  between  the  HULL  C  and  BAAL  runs  is  quite  remarkable, 
particularly  so  because  the  HULL  code  is  using  a  computational  grid 
designed  for  BAAL.  Normally,  the  HULL  code  uses  grids  in  which  the 
cell-to-cell  variation  in  size  is  kept  to  less  than  5  to  10  percent  in 
the  region  where  the  solution  is  desired,  or  in  any  region  through  which 
a  wave  must  pass  on  its  way  to  the  region  where  the  solution  is  desired. 
The  BAAL  grid  has  ccll-to-cell  variations  in  size  on  the  shelter  face 
of  25  percent  in  the  X  direction,  and  as  large  as  25  percent  for  most 
cells  in  the  Y  and  Z  directions,  with  much  larger  variations  in  the 
general  flow  field  away  from  the  shelter.  As  indicated  in  earlier 
work10  there  is  excellent  agreement  between  the  experimentally  measured 
average  overpressure  and  the  average  overpressure  computed  with  BAAL 
for  the  front  face  of  a  rectangular  parallelepiped,  using  a  similar  grid. 
This  agreement  gives  reason  for  confidence  in  the  validity  of  the  BAAL 
computation  for  the  shelter  front  face,  except  for  late  time  where  the 
average  overpressure  computed  with  BAAL  seems  high,  and  for  pressure 
anomalies  in  the  computation  at  the  edges  of  the  front  face. 

The  next  comparison  interest  in  Figure  7  is  that  of  HULL  C  with 
HULL  B.  The  input  data,  initial  flow  fields,  and  finite  difference 
grids  for  these  two  computations  are  identical  except  that  HULL  B  uses 
n  *  0.50  and  no  artificial  viscosity,  and  HUM.  C  uses  n  =  0.15  with 
artificial  viscosity  turned  on.  There  arc  significant  differences 
between  the  results  of  the  two  computations  early  in  the  diffraction 
phase.  Hie  HULL  R  computation  shows  an  overshoot  of  the  average  over¬ 
pressure  to  a  peak  value  of  88.5  kPa  at  1.0  ms,  12.7  percent  above  the 
peak  theoretical  value  of  78.5  KPa,  and  28.2  percent  above  the  peak 
average  overpressure  of  69.0  kPn  at  1.3  ms  in  the  HULL  C  computation. 

Hie  HULL  B  computation  also  shows  significant  oscillation  about  the 
results  of  the  HULL  C  computation  for  t  <  6  ms.  It  is  doubtful  that 
the  higher  CFl.  number,  n,  for  HULL  B  is  the  primary  cause  of  this 
oscillation;  it  is  more  likely  due  to  running  the  code  without  artifi¬ 
cial  viscosity,  combined  with  having  the  shock  started  as  a  disconti¬ 
nuity  at  the  shelter  front  face.  It  is  of  interest  to  mu  e  that  at 
1.6  ms  (the  HULL  R  datum  closest  in  time  to  the  first  datum  of  the  BAAL 
computation  at  1.571  ms)  the  IIUU.  B  computation  shows  an  average  over¬ 
pressure  of  76.4  kPa,  2.7  percent  less  than  the  theoretical  peak 
reflected  ow  .res sure  of  8.5  kPu. 


The  setup  of  the  HULL  A  computation  is  identical  with  that  of  the 
HULL  B  computation  except  that  the  HULL  A  computation  begins  with  the 
shock  located  seven  flow  field  cells  upstream  of  the  shelter  front  face 
As  can  be  seen  in  Figure  7,  the  movement  of  the  shock  through  those 
seven  flow  field  cells  results  in  a  significantly  different  average 
overpressure  history  on  the  shelter  front  face  for  HULL  A  compared 
with  HULL  B.  Signals  from  the  computational  shock  in  HULL  A,  spatially 
in  the  form  of  a  forward,  exponentially  decreasing  function,  arrive 
well  ahead  of  the  theoretical  shock  arrival  time.  These  early  signals 
are  then  reflected  back  continuously  into  the  main  portion  of  the  on¬ 
coming  computational  shock  wave,  ultimately  reducing  the  peak  over¬ 
pressure,  and  delaying  and  spreading  it  in  time.  The  peak  average 
overpressure  for  HULL  A  is  69.0  kPa,  12.1  percent  less  than  ti.e  theoret¬ 
ical  peak,  occurring  at  4.35  ms  after  the  theoretical  arrival  time  of 
the  shock  wave  at  the  front  face.  By  4.35  ms,  the  theoretical  location 
of  the  shock  wave  is  actually  1.68  m  downstream  of  the  front  face,  77.5 
percent  of  the  distance  to  the  plane  of  the  back  face.  Gentry  et  al.8 
show  the  same  qualitative  effect  in  comparing  the  results  of  a  BAAL 
computation  for  a  step  shock  started  at  the  front  face  of  a  rectangular 
parallelepiped  with  that  for  an  identical  shock  in  the  same  grid 
started  five  flow  field  cells  upstream  of  the  front  face. 

The  almost  exact  agreement  in  peak  average  overpressure  for  the 
shelter  front  face  between  HULL  C  and  A  is  most  likely  coincidental. 

The  effect  manifesting  itself  in  HULL  A  is  numerical  diffusion  caused 
by  using  a  CFL  number  less  than  1.0  (n  =  0.50)  and  having  the  shock 
move  through  seven  highly  nonuniform  flow  field  colls  prior  to  inter¬ 
acting  with  the  shelter.  The  X-direction  ccll-to-cell  variation  in 
size  between  the  seventh  and  sixth  cells  upstream  of  the  shelter  front 
face  is  45.0  percent.  The  early  time  effects  in  HULL  C  governing  the 
development  of  the  peak  average  overpressure  on  the  front  face  are  due 
almost  entirely  to  the  effects  of  the  modifications  via  artificial 
viscosity  to  tho  pressure  difference514*18  used  to  comnut e  the  accelera¬ 
tions  in  the  Lagrangian  phase  of  the  HULL  computation  .  The  HULL  A 
and  HULL  ft  computations  converge  with  one  another  by  20  ms,  as  should 
be  expected.  Tho  HULL  C  computation  appears  to  have  been  tending 
toward  converging  with  tho  other  HULL  computations  just  prior  to  its 
being  dostroyod  by  the  artificial  wave  discussed  earlier. 

The  solutions  obtained  by  using  the  model  outlined  in  the 
Technical  Manual  and  that  suggest  oil  by  Taylor  are  also  shown  in 
Figure  7.  All  six  average  overpressure  histories  suggest  that  the 
wave  interaction,  or  diffraction,  phase  on  the  front  face  lasts  for 
about  15  ms,  and  is  then  followed  by  the  drag  phase.  The  Technical 
Manual  solution  for  the  diffraction  phase  is  significantly  different 


0.  Rich tmycr  and  K.  V.  Morton.  "Difference  Method*  for  Initial 
Value  Problems !id  Rd>,  Tntarooienae  Publishers,  Tne.,  Rev  York 
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from  all  of  the  other  solutions,  except  that  it  is  nearly  identical  to 
Taylor's  solution  in  the  early  part  of  the  diffraction  phase.  It  agrees 
well  with  the  HULL  B  solution  throughout  the  drag  phase,  with  the  BAAL 
and  HULL  C  solutions  in  the  early  drag  phase,  and  with  the  HULL  A 
solution  shortly  after  the  beginning  of  the  drag  phase.  The  overpressure 
calculated  with  Taylor's  model,  although  it  is  somewhat  high,  shows  good 
agreement  with  the  hydrocode  solutions  throughout  the  diffraction  phase, 
except  for  regions  where  the  hydrocode  solutions  show  the  questionable 
behavior  discussed  earlier.  This  good  agreement,  coupled  with  the  ease 
with  which  the  model  may  be  applied,  makes  Taylor's  model  most  attractive 
for  estimating  front  face  diffraction  loading.  Taylor's  solution  for 
the  drag  phase  is  consistently  above  all  of  the  others,  except  for  the 
late-time  BAAL  solution. 

Appendix  D  contains  a  comparison  of  the  HULL  B  prediction  for 
average  overpressure  on  the  shelter  front  face  with  Ethridge's12 
empirical  correlation  function  and  a  gage  weighting  function  used  by 
Taylor1 . 

C.  Back  Face 


Figure  8  shows  a  comparison  of  the  average  overpressure  on  the 
back  face  of  the  S-280  shelter  as  computed  by  BAAL,  HULL  runs  A  and  B, 
the  Technical  Manual,  Taylor's1  suggested  improvement  to  the  Technical 
Manual  model,  and  a  proposed  model.  There  is  excellent  agreement 
among  the  three  hydrocode  computations.  The  results  for  the  two  existing 
semiempirical  models1*2  differ  from  those  for  the  hydrocodes,  particu¬ 
larly  the  model  from  the  Technical  Manual.  As  indicated  by  Gentry 
et  al.8,  the  average  overpressure  on  the  back  face  in  the  BAAL  computa¬ 
tion  for  the  rectangular  parallelepiped  agrees  quite  well  with  the 
experimentally  measured  average  overpressure.  This  led  to  the  conclusion 
that  the  average  overpressure  on  the  back  face  of  the  shelter  as 
computed  with  BAAL  may  also  bo  considered  to  be  an  accurate  estimate. 

This  conclusion  is  strongly  reinforced  by  the  agreement  shown  by  HULL 
runs  A  and  B.  The  three  hydrocode  computations  show  a  much  faster 
rise  time  required  to  roach  the  peak  average  overpressure,  and  that  the 
peak  is  higher  than  that  predicted  by  the  semiempirical  models.  The  . 
hydrocodes  also  indicate  a  local  pressure  peak  around  19  ms,  most 
likely  due  to  the  interaction  with  one  another  of  the  weakened  shocks 
breaking  over  the  top  and  around  the  side  faces. 

It  is  a  fairly  direct  matter  to  modify  the  existing  models  for 
average  overpressure  on  the  back  face  to  better  fit  the  hydroe ode 
results,  at  least  for  this  case.  As  before,  Equation  (3)  gives 
td  =  5.62  ms,  with  p*  »  0.0  for  t  <  tj.  It  is  proposed  that  the 

constant  a  in  Equation  (4)  bo  modified  so  that  n  =  2.  the  modified 
rise  time  is  then  (tb)risc  •  10.64  ms.  The  peak  average  overpressure 

on  the  back  face  may  be  computed  by  using  a  modified  version  of 
Equation  (a), 
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with  n  =  -  and  p*(t)^  =  (P2 -p^ 3 •  The  average  overpressure  on  the  back 

face,  p£,  varies  linearly  from  a  value  of  0.0  at  t  =  to  (Pb)max  at 

t  =  t ,  +  (t,  )  .  .  For  the  present  case,  (p?)  =  32.4  kPa  at  t  =  16.3  ms. 

d  b  rise  r  rb  max 


There  is  a  decompression  phase,  lasting  from 
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where  (p£)  varies  linearly  from  (Ppmax  to  the  value  (Ppdrag-  This 
value  is  found  by  substituting  n  =  0  into  Equation  (12),  so  that 
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Equation  (14)  gives  the  average  overpressure  on  the  back  face  for 
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a  step  shock  for  all  t  > 
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For  the  present  case, 


(i>.*)  =31.5  kPa  for  t  >  29.6  ms.  This  modified  model  is  shown  in 

vtbJdrag 

Figure  8.  It  fits  the  hydrocode  data  well,  but  its  generality,  like 
that  of  the  model  on  which  it  is  based,  may  be  questionable. 


D.  Top  Face 


Figure  4  shows  the  average  overpressure  on  the  top  face  of  the 
shelter  as  computed  by  BAAL,  HULL  runs  A  and  B,  the  Technical  Manual, 
and  a  proposed  model.  There  is  excellent  agreement  between  the  BAAL 
and  HULL  B  computations  throughout  the  range  of  simulated  time.  The 
HULL  A  computation  agrees  quite  well  with  the  other  two  hydrocode 
computations  for  0  <  t  <  7  ms,  and  for  t  >  16  ms.  For  7  <  t  <  16  ms, 
the  HULL  A  solution  shows  a  higher,  more  delayed  peak  compared  with 
the  other  computations.  This  is  most  likely  caused  by  the  delayed, 
extended,  and  loss  strong  reflection  that  the  HULL  A  computation  shows 
for  the  front  face.  The  HULL  A  computation  also  shows  some  top  face 
loading  for  t  <  0  ms,  indicating  the  arrival  of  tho  forward  portion  of 
the  computational  shock  ahead  of  the  theoretical  shock  arrival  time. 


Tho  solid  lino  in  Figure  4  represents  tho  average  overpressure  lor 
the  top  faco  as  computed  by  using  the  model  suggested  in  tho  Technical 
Manual.  Tho  agreement  with  the  hydrocodes  for  the  first  6  ms,  approxi¬ 
mately  equal  to  the  time  roquired  for  tho  theoretical  shock  wave  to 


travel  along  the  top  face,  is  excellent.  However,  the  model  under¬ 
predicts  the  peak  average  overpressure  computed  by  the  hydrocodes,  and 
shows  poor  agreement  for  t  >  6  ms.  As  with  the  models  for  the  back 
face,  the  Technical  Manual  model  for  the  top  and  side  faces  can  be 
modified  without  much  difficulty  to  give  excellent  agreement  with  the 
hydrocode  results. 


The  following  modified  version  of  Equation  (7)  is  proposed  for 
computing  the  average  overpressure  on  the  top  face, 
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where  A  is  a  constant  whose  value  depends  on  whether  the  maximum,  local 

minimum,  or  drag  phase  overpressure  is  to  be  computed.  After  a  delay 

time  t,  as  defined  in  Equation  (3),  the  maximum  average  overpressure 
d 

on  the  top  face,  (p*  )  ,  is  computed  by  setting  A  =»  1,  so  that 
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Equation  (15)  reduces  to 
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From  time  t  =  0  to  t  =  td,  the  average  overpressure  increases  linearly 
from  zero  to  (pJop)max<  Froir‘  time  t  '  td  t0  a  time 
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the  average  overpressure  on  tho  top  face  docroases  linearly  from 

(p*  )  to  a  value  (p*  .  which  is  computed  by  sotting  A  « 

top  max  1  top  pmin 

in  Equation  (15).  From  time  (tt0p)pro^n  to  a  time 
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the  average  overpressure  on  the  top  face  increases  linearly  from 

'PlonWn  t0  “  VOlU°  ®topW  “hlC"  ‘S  C0“PUt0<1  by  S'''"i''8  *  =  0'8 

in  Equation  (15);  after  this  time  tho  avorage  overpressure  remains  at 
(n*  1  Thus,  instead  of  dealing  with  several  different  equations 

V1  top'drug*  ,  .  , 

as  suggested  in  the  Technical  Manual7,  it  is  only  necessary  to  deal 

with  an  equation  of  the  form  of  Equation  (15),  with  A  =  1.0  for  the 
maximum,  A  =  0.6  for  the  minimum,  and  A  =  0.8  for  the  drag  phase 
average  overpressures.  For  the  the  present  cose,  the  model  proposed 
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above  predicts  p*o^  =0.0  kPa  at  t 


G.O  s,  (p*  )  =34.5  kPa  at 

’  top  max 


*d  ’  5’62  Bs’  ^top’pmin  “  26'7  kPa  a'  ^top’pmin  =  16'8  ”s'  a"d 

'PJop'drag  =  30'6  kPa  for  *  s  "top1  drag  ”  25'3  “•  The  histor>'  °f 
average  overpressure  for  the  top  face  as  predicted  by  the  modified  model 
proposed  above  is  shown  in  Figure  4. 


E.  Side  Face 


Figure  5  shows  the  average  overpressure  on  the  side  face  of  the 
shelter  as  computed  by  BAAL,  HULL  runs  A  and  B,  the  Technical  Manual, 
and  the  proposed  model.  The  proposed  model,  a  modification  of  the 
Technical  Manual  model,  is  also  to  be  used  for  both  the  top  and  side 
faces  of  the  shelter.  Because  of  the  geometry  of  the  shelter,  the  side 
face  is  nearly  equivalent  in  size  to  the  half  of  the  top  face  that  is 
actually  simulated  in  the  hydrocode  solutions  (the  top  face  is  cut  in 
half  by  a  symmetry  plane),  and  because  they  have  similar  orientations 
in  the  flow  field,  the  hydrocode  solutions  for  the  side  face  arc  nearly 
identical  with  those  for  the  top  face.  Therefore,  the  same  comments 
regarding  the  relative  performance  of  the  hydrocodes  and  the  existing 
models  for  the  top  face  loading  also  apply  to  the  side  face  loading. 

As  noted  for  the  modified  model  proposed  for  the  back  face  loading,  the 
generality  of  this  proposed  modified  model  may  also  be  limited. 


VII.  NET  LOADING 


A.  Definition 

An  important  point  stressed  by  Taylor1  is  that  the  method  outlined 
in  the  Technical  Manual  significantly  overestimates  the  net  loading  on 
a  3-D  target,  with  net  loading  defined  as  the  average  overpressure  on 
the  front  face  minus  that  on  the  back  face.  Taylor's  model  predicts  a 
lower  net  load  on  the  shelter  than  does  the  Technical  Manual  model . 


B.  Discussion 

It  is  appropriate  at  this  time  to  refer  back  to  the  BAAL  and 
HULL  B  plots  of  overpressure  versus  time  in  Figure  2  and  the  net  loading 
plots  shown  in  Figure  3.  There  are  several  reasons  why  only  the  BAA1. 
and  HULL  B  computations  arc  chosen  for  comparison.  The  HULL  (.  solution 
for  the  front  face  is  valid  only  up  to  19  ms,  as  discussed  earlier;  the 
solution  for  the  back  face  is  no  longer  available,  hut  in  view  of  the 
problems  associated  with  this  computation,  it  would  probably  be  of 
limited  value.  The  HULL  A  solution  shows  significant  deviation  ivom 
the  other  hydrocodc  solutions,  most  likely  due  to  the  diffusion  of  the 
computational  shock  after  its  having  been  started  seven  ce  Is  upstream 
of  the  shelter  front  face.  Both  the  BAAL  and  HULL  B  solutions  start 
with  the  computational  shock  at  the  front  face. 
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As  shown  in  Figure  2,  except  for  the  BAAL  computation  later  in  the 
drag  phase  (t  >  20  ms)  and  the  HULL  B  computation  near  t  =  0,  both 
hydrocode  solutions -are  consistently  below  the  front  face  average  over¬ 
pressure  predicted  by  either  Taylor’s  model  or  the  Technical  Manual 
model.  Conversely,  both  hydrocode  solutions  for  the  back  face  indicate 
that  the  back  face  is  loaded  even  more  rapidly  than  Taylor  suggests, 
peaking  and  remaining  at  average  overpressure  values  well  above  those 
predicted  by  either  of  the  two  models  for  all  of  the  time  represented 
here.  Viewed  together,  the  front  and  back  face  loadings  from  the 
hydrocode  solutions  show  a  greatly  reduced  net  load  on  the  shelter,  as 
may  be  seen  in  Figure  3.  The  curves  in  Figure  3  suggest  that  the 
diffraction  phase  for  the  entire  shelter  is  essentially  over  by  20  to 
25  ms,  or  approximately  four  crossing  times.  During  the  diffraction 
phase,  both  hydrocode  solutions  show  much  less  net  loading  than  that 
predicted  by  Taylor1,  which  is  in  turn  much  less  than  that  predicted  by 
using  the  model  in  the  Technical  Manual.  The  implication  is  that  the 
whole  body  response  of  a  3-D  target  during  the  diffraction  phase  (which 
can  be  relatively  long  compared  with  a  single  crossing  time  based  on  the 
dopth  of  the  structure)  may  be  significantly  less  than  the  simple  models 
predict . 

By  25  ms,  the  net  loadings  from  each  of  the  two  hydrocode  solutions 
shown  in  Figure  3  have  begun  to  diverge  from  one  another.  The  net 
loading  from  the  HULL  B  solution  increases  with  time,  hut  remains  less 
than  that  for  both  of  the  semicmpirical  model  solutions.  The  BAAL 
solution  eventually  becomes  greater  than  all  of  the  others;  this  may  he 
an  artifact  of  the  BAAL  computation,  but  that  is  yet  to  be  determined. 

The  data  in  Figure  5  suggest  that  the  drag  phase  for  the  structure  as  a 
whole  begins  at  20  to  25  as. 

Structural  damage  to  the  shelter  is  generally  done  during  the 
diffraction  phase*,  the  walls  are  driven  inward  and  can  damage  the 
equipment  inside,  reducing  or  ending  its  usefulness.  At  assumed  threat 
levels,  the  diffraction  phase  produces  relatively  little  whole-body 
motion  of  the  shelter;  this  motion  can  occur  in  the  drag  phase.  It  may 
be  that  the  shelter  and  the  equipment  inside  will  survive  the  diffraction 
phase  but  will  be  destroyed  by  having  the  shelter  roll  over  and  impact 
the  ground  in  the  drag  phase.  Thus,  understanding  what  the  shelter 
loading  is  during  the  drag  phase  is  essential.  Yet  it  is  during  the 
drag  phase  that  the  limited  computational  and  scmiempirical  data  that 
are  available  (such  as  are  shown  in  Figure  3)  disagree.  Additional 
experimental  and  computational  data  in  the  drag  phase  for  three- 
dimensional  structures  arc  needed  to  better  quantify  the  loading  during 
this  phase. 


VIII.  ROTATIONAL  MOMENT 


Neither  the  existing  semiempirical  models  nor  the  proposed  model 
provide  any  information  concerning  the  locations  of  the  centers-of- 
overpressure  on  the  surfaces  of  a  target  as  a  function  of  time;  wave 
interaction  models  could  provide  such  estimates.  The  TRUCK  code3,  which 
discretizes  surfaces  for  its  wave  interaction  computations,  does  provide 
this  capability.  To  the  extent  of  the  credibility  of  the  hydrocodes  and 
the  solutions  gained  by  their  use,  such  information  is  readily 
available1 0 . 

Figure  9  shows  the  history  of  the  Z  location  of  the  center-of- 
overpressure  for  the  front  face  of  the  shelter  as  computed  by  BAAL,  and 
by  HULL  runs  A  and  B.  The  three  hydrocode  solutions  show  remarkable 
agreement  with  one  another,  although  the  slight  differences  do  seem  to 
he  systematic.  The  center-of-overpressure  stays  quite  close  to  half 
way  up  the  front  face  for  nearly  all  of  the  simulated  time,  except  for 
a  slight  movement  downward  during  the  early  part  of  the  diffraction 
phase. 

Figure  10  shows  the  history  of  the  Z  location  of  the  center-of- 
overpressure  for  the  back  face  of  the  shelter  as  computed  by  the  same 
three  hydrocode  runs.  The  back  face  shows  considerable  variation  in 
the  location  of  the  center-of-overpressure  during  the  diffraction 
phase.  Prior  to  the  arrival  of  the  shock  wave,  the  cent er-of -over¬ 
pressure  is  arbitrarily  defined  as  being  at  one-half  of  the  height  of 
the  shelter;  during  that  time,  the  average  overpressure  is  zero  by 
definition.  The  HULL  A  plot  indicates  the  arrival  of  the  forward 
section  of  the  computational  shock  prior  to  the  theoretical  arrival 
time  of  t  =  5.62  ms  at  the  plane  of  the  hack  face. 

The  center-of-overpressure  histories  on  the  shelter  front  and  hack 
faces  take  on  added  significance  when  considered  in  conjunction  with 
the  average  overpressure  histories  on  the  respective  faces,  just  as  the 
average  loadings  on  the  front  and  back  faces  do  when  looked  at  together 
in  Figures  2  and  5.  The  semiempirical  models1 »?  predict  higher  loading 
on  the  front  face  than  do  the  hydrocodes.  For  lack  of  a  better  estimate, 
it  might  seem  reasonable  to  use  the  centroid  of  the  target  face  for  the 
location  of  the  center-of -overpressure  on  both  the  front  and  hack  faces 
when  using  the  semiempirical  models.  The  hydrocodes  indicate  that 
such  an  approximation  would  be  reasonable  for  the  front  face,  but  not 
so  for  the  diffraction  phase  on  the  back  face. 

During  the  early  part  of  the  diffraction  phase,  the  hydroeodcs 
predict  a  significantly  greater  average  overpressure  on  the  back  face 
than  do  the  semiempirical  models,  acting  considerably  above  the  mid¬ 
point  on  the  back  face.  During  the  remainder  of  the  time  simulated  by 
the  hydrocodes,  the  computed  louding  is  still  well  above  that  predicted 
by  the  semiempirical  models,  but  with  the  center -of -overpressure  near 
the  area  centroid.  Specifically,  an  average  of  the  BAAL  and  HULL  U 
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figure  9.  Height  of  centor-of -overpressure  on  front  face. 
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average  overpressure  on  the  hack  face  over  the  time  interval 
6.6  <  t  <  IS. 8  ms  is  1S1  percent  greater  than  that  predicted  by  the 
Technical  Manual,  and  S7  percent  greater  than  that  predicted  by  Taylor's 
model.  For  the  time  interval  17.2  <  t  <  33.4  ms,  the  two  hydrocode 
solutions  average  30  percent  greater  than  the  Technical  Manual  and  8.7 
percent  greater  than  Taylor's  model. 


The  histories  of  the  net  rotational  moment  due  to  overpressure  on 
the  shelter  computed  by  using  BAAL,  HULL  8,  and  each  of  the  two  semi- 
empirical  models  are  shown  in  Figure  11.  The  rotational  moments  are 
computed  about  a  line  in  the  plane  of  the  bottom  boundary  and  perpendic¬ 
ular  to  the  side  face.  As  was  the  case  for  the  net  loading  curves  in 
Figure  3,  Taylor's  model1  yields  rotational  moments  larger  than  those 
for  the  Technical  Manual  very  early  in  the  diffraction  phase  and  also 
in  the  drag  phase,  but  significantly  smaller  values  throughout  the 
majority  of  the  diffraction  phase.  The  HULL  B  computation  shows  the 
oscillation  in  the  solutions  in  the  early  diffraction  phase  as  discussed 
earlier.  For  4  <  t  <  20  ms  the  HULL  B  and  BAAL  solutions  show  excellent 
agreement  with  one  another,  predicting  values  well  below  those  for  the 
simple  models.  During  that  time  interval,  the  hydrocodes  average  47 
percent  less  than  that  predicted  by  the  Technical  Manual,  and  '>1 
percent  less  than  that  predicted  by  Taylor's  mtxlel .  For  t  >  20  ms,  the 
BAAL  and  HULL  B  solutions  diverge  significantly  from  one  another.  By 
33.1  ms,  the  HULL  B  computation  has  reached  a  net  rotational  moment  due 
to  overpressure  of  55.0  kN-m,  20.3  percent  less  than  the  Technical 
Manual  model's  value  of  69.0  KN-m,  and  33.5  percent  less  than  the  value 
of  82.7  kN-m  obtained  by  using  Taylor's  model .  By  33.6  ms,  the  RAAf 
computation  has  reached  a  net  rotational  moment  due  to  overpressure  of 
97.8  kN-m,  41.7  percent  greater  than  the  value  from  the  Technical 
Manual  model,  19.5  percent  greater  than  the  value  from  Taylor's  model, 
and  77. s  percent  greater  than  the  value  from  the  MUM.  B  computation  at 
33.1  ms. 


The  reason  for  this  divergence  is  two-fold.  The  BAM  computation 
shows  much  higher  late-time  overpressure  on  the  front  face  of  the 
shelter  than  does  the  MUM,  B  computation,  as  shown  in  Figure  7.  with  a 
center -of -overpressure  consistently  above  that  for  the  HUM.  B  computa¬ 
tion,  as  shown  in  Figure  9.  Conversely,  the  BAM  computation  shows 
only  somewhat  higher  late-time  overpressure  on  the  hack  face  of  the 
shelter  than  does  the  HUM,  B  computation,  as  shown  in  Figure  s* 
a  center -of -overpressure  consistently  below  that  for  the  HUM  B  compu¬ 
tation,  as  shown  in  Figure  JO.  At  least  part  of  the  reason  for  this 
difference  in  behavior  between  the  two  hydroeodc  solutions  is  the 
apparent  anomalous  behavior  in  the  BAAL  solution  discussed  earlier 
More  information  on  the  grid  and  the  anomalous  behavior  is  presented 
in  Appendices  C  and  E. 


If  the  hvdrocodc  computations  are  correct,  as  they  appear  to  he, 
then  the  net  rotational  moment  on  a  3-D  structure  may  he  much  less  than 
would  be  indicated  by  the  existing  seta i empirical  models.  This  uuiu^n 
loss  vulnerability  to  overturning  than  may  have  previously  been  assumed 
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IX.  CONCLUSION 


In  general,  the  hydrocodes  show  good  agreement  with  one  another  in 
predicting  overpressure  averaged  over  each  of  the  given  faces  of  the 
S-280  Electrical  Equipment  Shelter.  They  are  also  useful  as  comparative 
tools  in  evaluating  the  capability  of  the  semi  empirical  models  to 
predict  such  loading.  Additionally,  the  hydrocodes  provide  a  great  deal 
of  information  for  which  more  easily  used  predictive  models  (similar  to 
the  semi empirical  models  discussed  here)  are  not  available  and  which  is 
not  readily  gathered  experimentally.  An  example  of  this  is  the  use  of 
the  data  from  the  hydrocodes  to  make  accurate  estimates  of  the  location 
vs  time  of  the  cent er-of -overpressure  for  use  in  computing  rotational 
moments.  Also,  the  analysis  of  the  various  hydrocode  results  discussed 
here  has  allowed  some  quantification  of  the  effects  of  such  items  as 
the  value  of  the  CFL  number,  numerical  diffusion  caused  by  the  finite 
difference  grid,  and  artificial  viscosity. 


The  model  outlined  in  the  Technical  Manual  does  not  appear  to  be 
suitable  for  predicting  either  the  shock  diffraction  or  the  drag  phase 
loading  on  a  three-dimensional  structure  such  as  the  S-280  Electrical 
Equipment  Shelter.  The  front  face  loading  prediction  for  the  diffraction 
phase  seems  to  be  too  high;  the  prediction  for  the  drag  phase  appears 
to  be  reasonable.  The  prediction  for  the  back  face  loading  is  too  low, 
and  is  made  worse  by  predicting  too  long  a  rise  time  for  overpressure 
to  go  from  zero  to  the  predicted  drag  phase  value.  The  models  for 
predicting  the  loading  on  the  top  and  side  faces  are  fairly  good  at 
predicting  the  initial  increase  of  average  overpressure  with  time,  but 
they  underestimate  its  peak  value;  thereafter,  the  models  uppear  to  be 
significantly  in  error. 


Taylor's  model1  for  estimating  the  average  overpressure  on  the 
front  of  a  thrcc-dimonsionol  structure  is  a  significant  improvement 
over  the  model  suggested  in  the  Technical  Manual,  particularly  during 
the  diffraction  phase.  However,  predictions  of  average  overpressure 
in  the  drag  phase  using  Taylor's  model  appear  to  he  somewhat  high 
compared  with  the  hydrocode  predictions.  The  change  suggested  by 
Taylor  to  the  Technical  Manual  model  for  the  average  overpressure  on 
the  back  face  is  an  improvement,  but  it  too  underprediets  the  loading 
for  all  time,  missing  the  peak  overpressure  and  overestimating  the 
timo  required  to  reach  the  predicted  peak  overpressure.  Taylor  did 
not  address  the  top  and  side  face  loading. 


The  empirical  models  suggested  by  Ethridge1'  provide  a  quick  and 
fairly  accurate  means  of  estimating  the  average  ovorpressure  as  a 
function  of  time  on  both  the  front  and  back  faces  of  a  rectangular 
parallelepiped.  The  modifications  to  the  models  in  the  Technical 
Manual  for  the  hack,  top,  and  side  faces  of  a  rectangular  parallelepiped 
proposed  in  this  report  also  provide  a  quick,  accurate  means  of  esti¬ 
mating  average  overpressure,  with  the  added  feature  of  providing  some 
of  the  detailed  loading  variations  observed  in  the  hydroeodc  computa¬ 
tions. 


42 


The  semiempirical  models  do  not  provide  a  means  for  computing  the 
time-history  of  the  center-of-overpressure  on  the  various  faces  of  th< 
structure.  As  a  result,  it  is  not  possible  to  compute  rotational 
moments  unless  assumptions  arc  made  concerning  the  location  of  the 
center-of-overpressure.  Alternatively,  center-of -overpressure  models 
could  he  developed.  The  hydrocodes  do  provide  this  information. 

The  hydrocodes  and  the  semiempirical  models  all  provide  conflicting 
information  during  the  drag  phase.  This  late-time  loading  is  important 
because  it  is  during  this  time  that  overturning  either  will  or  will  not 
take  place.  Unfortunately,  relatively  little  experimental  data  have 
been  gathered  for  late-time  loading.  This  problem  needs  to  be  resolved 
both  experimentally  and  computationally.  Comparison  of  late-time 
loading  obtained  via  experiment  with  that  from  inviseid  codes  such  ns 
HULL  will  also  be  of  value  in  determining  whether  or  not  viscous  effects 
should  be  modeled. 
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Al'PliNDIX  A 


SOMITION  !;OR  SHOCK  TUBfi  WITH  Cl.OSI-D  »:Nl> 


i 


The  following  is  a  summary  of  the  flow  conditions  , or  < he  MS  kPa 
(5.0  psi)  step  shock  described  in  the  section  "Reference  Problem"  and 
used  in  the  HULL  and  BAAL  computations.  The  shock  is  assumed  to  have 
been  produced  in  a  shock  tube  having  a  high-pressure  reservoir  and  a 
low-pressure  reservoir  initially  in  thermal  equilibrium  with  one 
another,  with  the  fluid  in  each  reservoir  at  rest.  The  high-pressure 
reservoir  corresponds  to  region  3  in  Figure  1,  and  the  low-pressure 
reservoir  to  region  1.  The  shock  tube  is  assumed  to  be  closed  at  both 
ends. 

At  some  time,  t',  after  the  diaphragm  separating  the  high-  and 
low-pressure  reservoirs  has  been  removed,  a  right -travel ing  shock  wave 
has  developed  and  is  moving  into  region  1,  producing  shocked  region  2. 
Simultaneously,  a  left-traveling  expansion  fan  has  developed  and  is 
moving  into  region  3,  producing  the  expanded,  accelerated  region  4. 
Regions  2  and  4  are  adjacent  to  one  another.  They  have  matching  pressure 
and  normal  velocity  components  across  their  common  boundary,  but  will, 
in  general,  have  different  densities  and  specific  internal  energies. 
Hence,  this  boundary  is  a  discontinuity,  called  a  contact  discontinuity 
or  contact  surface.  At  some  time,  t",  the  shock  reflects  from  the 
closed  right  end  of  the  shock  tube,  creating  a  left-traveling  reflected 
shock  wave  with  region  5  behind  it. 

For  the  reference  problems  T^  =  288.2  °K,  pj  -  101.3  kPa,  and 

pj  =  1.225  kg/nf*.  The  gas  in  both  regions  1  and  3  is  assumed  to  be  air, 

and  is  further  assumed  to  be  a  polytropic  gas  with  >  =  1.4.  The  shock 
has  an  overpressure  of  34,5  kPa.  Table  A.l  contains  a  fairly  complete 
description  of  the  five  regions,  using  numerical  values  taken  to 
slightly  higher  accuracy  than  the  nominal  values  used  in  the  text. 

The  computational  grid  was  set  up  using  the  conditions  in  region  1 
for  the  undisturbed  ambient  air,  and  those  in  region  2  for  the  shocked 
air.  The  conditions  in  region  5  represent  the  region  behind  the 
reflected  shock  prior  to  the  arrival  of  any  expansion  waves 


APPENDIX  B 

DETAILED  DISCUSSION  OF  SEMIEMPIRICAL  MODELS 
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A.  Front  Face 


The  Technical  Manual's  method  for  computing  the  loading  on  the 
front  face  of  a  rectangular  parallelepiped  (such  as  the  shelter)  is 
based  on  an  assumed  clearing  time  t  ,  where 


t  = 
c 


(B.l) 


Here,  c^  is  the  sound  speed  on  the  front  face  of  the  shelter  after  the 

shock  has  reflected  from  it.  The  value  to  be  used  for  h  is  the  smaller 
of  either  the  height  of  the  shelter  or  one  half  the  width.  The  clearing 
time  t  is  "...  the  time  required  to  clear  the  front  wall  of  reflection 

effects"2.  During  this  time,  the  average  overpressure  on  the  front  face 
decreases  from  the  peak  reflected  overpressure  to  a  value,  in  the  case 
of  an  exponentially  decaying  wave,  computed  by2 

Pf  «  p*(t)f  +  0.85  j  p2  u22.  (B.2) 


In  F.quation  (B.2),  p2  is  the  density  behind  the  incident  shock,  u9  is 

the  particle  velocity  behind  the  incident  shock,  and  p*(t)f  is  the  time- 

dependent  value  of  the  incident  shock  overpressure  at  the  plane  of  the 
front  face, 


p*(t){  =  » -  (t .  tdyti)os 


(B.3) 


Here,  pj  is  the  absolute  pressure  ahead  of  the  shock  and  p2  is  the 


absolute  pressuro  behind  the  shock.  Also,  t  is  the  time  measured  after 
shock  arrival  at  the  front  face;  t(1  is  the  time  required  to  travel  to 

the  surface  under  consideration,  which  in  this  case  is  the  front 
surfaco,  so  t^  =  0.0;  tpQs  is  the  duration  of  the  positive  phase.  For 

a  stop  shock,  Equation  (B.3)  rcducos  to 


P*(t)f  B  Pf  °  (P2  -  Pp* 


(B.4) 


so  Equation  (B.2)  becomes 

Pf  ° 


1  2 

p2  -  p,  ♦  0.85  jp,  u2  . 


(B.5) 


For  the  reference  shock  (Sec  Section  II  and  Appendix  A.),  the  peak 
reflected  overpressure  on  the  front  face  of  a  rectangular  parallelepi¬ 
ped  is  (p5  -  pj). 
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Taylor1  suggested  an  alternate  method  to  estimate  the  loading  on 
the  front  face  of  an  obstacle,  based  on  the  number  of  rarefaction  wave 
crossings.  For  a  3-D  object  such  as  the  shelter,  there  will  be  a 
succession  of  rarefaction  waves  originating  from  the  top  of  the  front 
face  and  reflecting  from  the  bottom  boundary,  plus  a  succession  of 
rarefaction  waves  originating  from  the  side  face  and  reflecting  from 
the  symmetry  plane.  The  reflected  waves  will  also  be  expansion  waves. 

In  an  actual  occurrence  of  loading  through  shock  diffraction,  there  is 
considerable  interaction  between  the  crossing  expansion  waves.  Because 
Taylor  intended  that  the  method  be  a  simple  approximation,  the  waves 
are  assumed  not  to  interfere  with  one  another.  As  a  further  simplifi¬ 
cation,  the  rarefaction  wave  speeds  are  assumed  to  be  equal  to  the 
sound  speed  in  the  reflection  region  immediately  after  the  incident  shock 
reflects  from  the  front  face. 

Using  Taylor's  method,  the  required  crossing  time  for  an  expansion 
wave  running  from  the  top  to  the-  bottom  of  the  front  laco  is  computed, 
as  is  that  for  an  expansion  wave  t'unning  from  the  side  edge  to  the 
symmetry  plane.  For  each  expansion  wave  crossing,  the  average  over¬ 
pressure  on  the  front  face  is  then  computed  as  a  percentage  of  the 
initial  reflected  overpressure  on  the  front  face  (p<.  -  pj)»  as  shown  in 
Table  1  of  this  report.  "  . 

R,  Back  Face 

The  time  required  for  the  incident  shock  wave  to  arrive  at  the 
plane  of  the  back  wall  is 

-  §7  •  (B-« 

whore  D  is  the  shelter  depth,  and  Wj  is  the  velocity  of  the  incident 

shock.  The  rise  time*  required  for  the  pressure  on  the  back  face  to 
go  from  ambient  at  time  to  its  peak  average  value  is 


(B.7) 


where  a  ®  4.0,  Cj  is  the  sound  speed  ahead  of  the  shock,  and  h  is 
defined  as  before.  The  peak  average  overpressure*  on  the  back  face  is 

(hwi|'(,h(!)[i  * "  •  8,c‘6]  •  '"•R 

where  p*(t)b  is  the  incident  shock  wave  overpressure  at  the  time  it 
arrives  at  the  plane  of  the  back  face  (See  Equation  (B.S).),  and 


(B.9) 


Pi 


Here  [p*(t)^]t_o  is  the  incident  shock  wave  overpressure  at  the  time  it 
arrived  at  the  plane  of  the  front  face,  and  p^  is  the  absolute  pressure 
ahead  of  the  shock.  For  the  step  shock  considered  here, 

[P*(t)flt_0  =  P*(t)b  -  (p2  -  Pi)  (B-10) 

The  Technical  Manual  model  modified  for  a  step  shock  implies  that  the 
average  overpressure  on  the  back  face  would  remain  at  (ppmax  indefi¬ 
nitely.  * 

Taylor's  data  indicate  that  the  constant  a  used  in  Equation  (B.7) 
to  compute  the  rise  time  should  be  changed  to  2.5  from  4.0.  The  rise 
time  is  measured  from  the  arrival  time,  t^,  of  the  shock  at  the  plane 

of  the  back  face  to_the  time  at  which  the  average  overpressure  on  the 
back  face  reaches  (ppmax*  This  gives  a  more  rapid  rise  time  on  the 

back  face  than  that  computed  by  using  the  model  in  the  Technical 
Manual . 


C.  Top  and  Side  Faces 


Taylor  does  not  suggest  an  improved  method  for  computing  the 
average  overpressure  on  the  top  and  side  faces  of  an  object  such  as  the 
shelter.  The  method  outlined  "in  the  Technical  Manual  is  presented  here 
for  later  comparison  with  the  results  computed  using  the  hydrodynamic 
computer  codes.  From  time  zero,  when  the  incident  shock  wave  reaches 
the  plane  of  the  front  face  of  the  target,  to  the  time  t^  at  which  the 

shock  wave  reaches  the  plane  of  the  back  face,  the  average  overpressure 
on  the  top  face  varies  linearly5  from  zero  to 


0.9  ♦  0.1 


(tun 


for  the  step  shock  considered  here.  To  account  for  vortex  growth, 
shedding,  and  subsequent  travel  down  the  top  face,  the  model  predicts 
a  local  minimum  in  the  average  overpressure  at  the  time 


t  *  SI1 

pain  tT. 


(0.121 


At  that  time,  the  average  overpressure  on  the  top  face  is  computed-  as 
either 


S4 


Because  the  HULL  and  BAAL  codes  use  different  conventions  for 
naming  directions  and  indices,  the  BAAL  results10  and  the  flow  field 
description  are  reproduced  in  part  using  the  HULL  hydrocode  conventions. 
Figure  6  of  this  report  shows  a  3-D  view  of  the  shelter  in  the  computa¬ 
tional  flow  field.  Figure  C.l  shows  the  projection  of  the  computational 
grid  on  the  front  face.  Figure  C.2  the  projection  on  the  top  face,  and 
Figure  C.3  the  projection  on  the  side  face.  The  indicated  pressure 
anomalies  on  these  faces  in  the  BAAL  computation  are  discussed  in 
Appendix  E. 

The  computational  grid  used  for  the  HULL  code  is  identical  with 
that  used  in  the  BAAL  computation,  except  that  the  HULL  grid  has  one 
additional  plane  of  cells  in  the  X  direction,  four  additional  planes  of 
cells  in  the  Y  direction,  and  three  additional  planes  of  cells  in  the 

Z  direction,  all  added  at  the  upper  extreme  of  distance  in  the  respective 

directions.  Table  C.l  is  a  tabulated  list  of  ceU  vertex  locations 
versus  computational  grid  indices.  The  single  plane  of  cells  was  added 

to  the  HULL  grid  in  the  X  direction  to  satisfy  a  code  requirement  that 

the  number  of  cells  in  the  X  direction  be  odd.  The  extra  planes  of 
cells  were  added  in  the  Y  and  Z  directions  for  the  HULL  grid  to  allow 
the  code  to  process  a  longer  simulated  flow  time  prior  to  the  possible 
arrival  at  the  shelter  of  artificial  signals  reflecting  from  che 
transmissive  boundaries  of  the  computational  flow  field. 
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Figure  C.l.  Computational  grid  on  front  face. 

(Identical  back  face,  including  pressure 
anomalies.) 
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Table  C.l. 

Cell  Vertices 

,  HULL  and 

BAAL  Grids 

I 

X  (m) 

J 

Y  (m) 

K 

Z  On) 

Origin 

0.0000 

Origin 

0.0000 

Origin 

0.0000 

1 

1.9622 

1 

0.3376 

1 

0.3933 

2 

3.3154 

2 

0.6752 

2 

0.7865 

3 

4.2487 

3 

1.0128 

3 

1.1798 

4 

4.8923 

4 

1.2828 

4 

1.4944 

5 

5.3520 

5 

1.4989 

5 

1.7461 

6 

5.6804 

6 

1.6717 

6 

1.9474 

7 

5.9233 

7 

1.8100 

7 

2.1085 

8 

6.1176 

8 

1.9483 

8 

2.2696 

9 

6.2731 

9 

2.1332 

9 

2.4850 

10 

6.4285 

10 

2.3805 

10 

2.7732 

11 

6.6228 

11 

2.7383 

11 

3.1900 

12 

6.8657 

12 

3.2744 

12 

3.8145 

13 

7.1693 

13 

4.1965 

13 

4.7499 

14 

7.5488 

14 

5.6671 

14 

6.2012 

15 

7.8524 

15 

8.0129 

15 

8.4529 

16 

8,0953 

16* 

10.3587 

16* 

10.7047 

17 

8.2896 

17* 

12.7044 

17* 

12.9564 

18 

8,4451 

18* 

15.0502 

18* 

15.2081 

19 

8.6005 

19* 

17.3959 

20 

8.7948 

21 

9.0377 

22 

9.3661 

23 

9.8258 

24 

10.4694 

25 

11.4027 

26 

12.7559 

27 

14.7181 

28 

16.6803 

29 

18.6425 

30 

20.6047 

31* 

22.5668 

S 


Added  cell*  HULL  grid  only 
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APPENDIX  D 

OTHER  LOADING  FUNCTIONS 


Figure  D.l  shows  the  gage  positions  that  Taylor1  used  in 
estimating  the  loading  for  the  front  and  back  faces  of  a  three- 
dimensional  rectangular  parallelepiped  tested  in  the  BRL  shock  tube. 

The  experimentally  measured  overpressure  histories  for  these  points  were 
used  as  a  basis  for  evaluating  the  BAAL  computation  by  Gentry  et  al8. 
Point  A  is  located  at  the  centroid  of  triangle  FGO,  point  C  at  the 
centroid  of  triangle  EOH,  and  point  D  at  the  centroid  of  triangle  EFO. 
Taylor's  function  for  estimating  the  average  overpressure  on  either  the 
front  or  back  face  is 


P*  = 


2pA  +  PS  +  PD 
- 3 - 


(D.l) 


where  p*  is  the  overpressure  at  point  A,  p*  is  the  overpressure  at 
point  C,  and  p*  is  the  overpressure  at  point  D. 


Ethridge1 2  has  proposed  empirical  correlation  functions  for  the 
front  and  rear  surfaces  of  a  rectangular  parallelepiped.  For  the 
average  overpressure  on  the  front  surface  of  a  S-280  shelter  struck  by 
the  reference  34.5  kPa  step  shock,  Ethridge's  function  reduces  to 


pf 


38.51 


*  (40.02)  e<--!28S><374-4  0(I'412* 


.1667  e' 


l) 


• (D. 2) 


Figure  D.2  shows  a  comparison  of  the  average  overpressure  on  the 
front  face  computed  with  Ethridge's  empirical  function,  HULL  B,  and 
Taylor’s  weighted  function  shown  in  Equation  (D.l).  The  overpressure 
histories  for  points  A,  C,  and  D  were  computed  using  spatially  inter¬ 
polated  values  for  those  points  in  the  HULL  B  computation.  Figure  D.2 
indicates  that  Taylor's  woighted  function  for  those  three  points  gives 
an  excellent  approximation  to  the  average  overpressure  for  the  front 
face.  Ethridge's  empirical  function  also  shows  good  agreement  with  the 
computed  solution,  especially  at  late  time.  At  early  time,  F.thridgo's 
function  is  probably  more  reliable  than  the  HULL  B  solution  because  it 
is  based  on  the  theoretical  reflected  value  at  time  zero  and  is  not 
subject  to  the  kind  of  numerical  oscillation  exhibited  by  the  HULL 
solution.  However,  it  does  not  show  the  apparent  over-relief  of 
pressure  found  both  experimentally  and  computationally  near  10  ms. 

Figure  D.3  shows  the  average  overpressure  on  the  back  face  of  the 
shelter  computed  with  Ethridge's  empirical  function,  Taylor's  weighted 
function,  and  HULL  B.  Ethridge's  empirical  function  for  the  average 
overpressure  on  the  back  face  of  the  shelter  struck  by  the  reference 
shock  reduces  to 


(.9300)  1 


* 


c-161.2  t  (1  ♦  5656  t2)l 


(D.3) 
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Taylor's  weighting  function,  Equation  (D.l),  is  the  same  for  the  back 
face  as  for  the  front  face.  The  overpressures  for  points  A,  C,  and  D 
on  the  back  face  were  computed  by  performing  linear  interpolations  in 
the  HULL  B  solution  in  the1  same  manner  as  was  done  for  the  corresponding 
points  on  the  front  face.  Ethridge's  empirical  function  shows  good 
general  agreement  with  the  HULL  B  computation,  but  does  not  show  the 
peak  and  local  minimum  suggested  in  Figure  8  of  this  report  and  by  the 
proposed  model.  It  also  seems  to  load  the  back  face  more  rapidly.  The 
weighted  function  of  Equation  (D.l)  seems  to  be  a  reasonable  approxima¬ 
tion  to  the  full  solution  except  that  it  overpredicts  the  loading  near 
the  peak.  This  overpr edict ion  is  especially  interesting  because  it  is 
of  essentially  the  same  nature  as  the  overpr edict ion  of  Taylor's 
weighted  function1,  as  shown  by  Gentry  et  al8.  Figure  14  of  their 
report,  reproduced  here  as  Figure  D.4,  compares  the  experimentally 
measured  overpressures,  the  complete  BAAL  solution,  and  the  interpolated 
values  for  analogous  points  A,  C,  and  D  on  a  rectangular  parallelepiped 
tested  in  the  BRL  shock  tube.  The  solid  line  represents  the  weighted 
average  for  the  experimentally  measured  values  at  points  A,  C,  and  D. 

The  dashed  line  with  the  x's  represents  the  weighted  average  of  the 
overpressures  at  points  A,  C,  and  D  found  by  performing  linear  inter¬ 
polation  using  a  BAAL  hydrocode  computation.  The  dashed  line  with  the 
o's  represents  the  average  overpressure  on  the  back  face  using  the  BAAL 
computation  for  the  entire  face. 
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E  F 


Taylor's  gage  positioning  for  front  and 
back  face  loading  on  a  rectangular 
parallelepiped. 
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face  average  overpressure 
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Figure  D.4.  Average  overpressure  on  the  hack  face  of  a 
rectangular  parallelepiped  being  struck  by 
a  S.O  psi  step  shock.  (Reproduced  from 
Gentry  et  al0.) 


APPENDIX  E 

DETAILED  COMPARISON  OF  HYDROCODE  RESULTS 
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A.  Front  Face 


As  may  be  seen  by  comparing  Figures  C.l,  C.2,  and  C.3  in  Appendix 
C  with  Figures  2a,  2b,  and  2c  in  the  earlier  study10,  additional 
apparent  anomalies  have  been  identified  in  the  BAAL  computation  during 
the  current  analysis.  The  anomalies  on  the  front  face  of  the  shelter 
are  manifested  as  excessively  high  pressures  in  the  topmost  row  of  flow 
field  cells  (indices  I  =  9,  .1  =  1  through  7,  K  =  7,  or  (9,1-7,71)  and 
the  right-most  column  of  cells  (9, 7, 1-7).  Figure  F..1  shows  over¬ 
pressure  histories  for  cell  (9,7,4),  a  flow  field  cell  in  the  right¬ 
most  column  of  cells  on  the  shelter's  front  face,  four  cells  up  from 
the  bottom.  Except  for  very  early  in  the  diffraction  phase,  the  BAAl. 
overpressures  for  the  cells  in  that  column,  as  cell  (9,7,4)  typifies, 
are  in  general  significantly  larger  than  those  for  the  two  HULL  compu¬ 
tations.  Correspondingly,  Figure  E.2  shows  an  even  more  severe  effect 
in  the  overpressure  history  for  ceil  (9,7,7),  the  flow  field  cell  on 
the  upper  right  comer  of  the  front  face.  P.xcept  in  those  areas  where 
apparent  anomalous  behavior  is  exhibited  in  the  BAAl.  solution  for  the 
front  face,  the  relative  agreement  in  the  BAAL  and  HULL  computations 
is  more  accurately  characterized  by  the  overpressure  history  of  cell 
(9,3,3),  a  cell  located  near  the  center  of  the  one-half  shelter  front 
face,  as  shown  in  Figure  E.3.  Here,  the  BAAL  computation  shows  the 
more  typical  increase  in  overpressure  above  that  for  the  HULL  computa¬ 
tions  for  late  time  (t  >  20  ms). 

B .  Back  Face 

The  pressure  anomalies  in  the  BAAL  computation  for  the  hack  face 
of  the  shelter  are  in  essentially  the  opposite  sense  from  those  on  the 
front  face.  The  overpressures  for  the  top  row  of  flow  field  cells 
(19.1-7,7)  and  the  outermost  column  of  flow  field  cells  (19,7,1-7) 
are  significantly  less  than  those  computed  by  using  KUtl..  Figure  F..4 
shows  the  overpressure  history  of  cell  (19,3,7),  one  of  the  flow  field 
cells  along  the  top  row  of  cells  on  the  shelter's  hack  face.  It  shows 
the  typical  late-time  overpressure  disagreement  between  the  BAAl- 
computation  and  the  HULL  computations  for  that  row  of  flow  field  cells. 
Figure  K.S  shows  a  similar,  typical  disagreement  between  tbe  BAAL  and 
HULL  computations  for  a  flow  field  cell  in  the  outermost  column  of 
cells  on  the  back  face.  Figure  U.b  shows  tbe  overpressure  history  of 
flow  field  cell  (19,3,3),  situated  approximately  on  the  center  of 
the  one -ha  If  back  face.  The  relative  performance  of  the  BAAL  and  HULL 
computations  shown  in  Figure  F,.b  typifies  that  for  th  flow  field  cells 
on  the  back  face,  but  not  at  the  outer  edges  of  tbe  h-sck  face. 

6.  Top  Face 

Apparent  pressure  anomalies  in  tbe  BAAl.  computation  also  exist  on 
the  top  and  side  faces,  and  are  quite  similar  to  those  on  the  front 
and  back  faces.  Figure  -H.7  shows  the  overpressure  history  of  flow 
field  cell  (10,3,8),  the  third  cell  in  the  Y  direction  from  tbe 


symmetry  plane  in  the  first  row  of  flow  field  cells  along  the  top  face. 
Here,  the  overpressure  as  computed  by  BAAL  is  consistently  much  less 
than  that  computed  by  either  of  the  two  HULL  computations.  The  over¬ 
pressure  history  of  flow  field  cell  (9,5,7),  the  top  flow  fir'd  cell 
in  the  third  column  of  cells  up  the  front  face  (not  shown  in  this 
report),  shows  relative  values  of  overpressure  between  the  BAAL  and  HULL 
-..Computations  similar  to  those  shown  in  Figures  E.4  and  F..5.  This  is 
part’  •0-f.a  consistent  pattern  of  the  main  differences  in  the  overpressure 
histories  in-'-the  HULL  and  BAAL  computations  for  flow  field  cells  just 
before  and  just  after-  .a  90-degree  expansion  corner.  The  BAAL  computation 
shows  consistently  higher  values  of  overpressure  than  do  the  HULL  compu¬ 
tations  in  the  flow  field  cells' on  .the  shelter  surface  just  upstream  of 
the  90 -degree  expansion  corner,  and  consistently  lower  values  on  the 
shelter  surface  just  after  the  expansion.  Figure.  E. 8  shows  the  over¬ 
pressure  history  of  flow  field  cell  (18,3,8),  the  third  cell  in  the 
Y  direction  from  the  symmetry  plane  in  the  last  row  of  flow  field  cells 
along  the  top  face.  Here  the  overpressure  from  the  BAAL  computation'  is.... 
consistently  greater  than  from  the  HULL  computations.  Figures  F..1  and 
E.2  show  similar  results  for  analogous  front  face  cells.  Figure  E.9 
shows  the  overpressure  history  of  flow  field  cell  (14,7,8),  the  outer¬ 
most  cell  in  the  fifth  row  of  flow  field  cells  down  the  top  face.  This 
29-degree  corner  of  the  shelter  is  along  a  line  parallel  rather  than 
perpendicular  to  the  flow;  hence  there  are  no  large  pressure  or  velocity 
gradients  around  this  cell  as  there  would  be  for  a  similar  corner  which 
is  oriented  perpendicular  to  the  flow.  As  may  be  seen  in  Figure  f.9, 
the  agreement  between  the  HULL  and  BAAF  cedes  is  generally  good,  except 
for  the  delayed,  enhanced  peak  for  HULL  A;  the  cause  of  this  discrepancy 
was  discussed  earlier  in  this  report.  Figure  F..19  for  flow  field  cell 
(14,3,8),  located  near  the  center  of  the  half  of  the  top  face,  shows 
the  typical  good  agreement  between  the  HULL  and  BAAL  computations  for 
flow  field  cells  on  the  top  face  but  not  immediately  adjacent  to  the 
90-degree  expansion  comer  regions. 

D.  Side  Face 


Figure  fi.ll  shows  the  overpressure  history  of  flow  field  cell 
(10,8,4),  the  fourth  cell  up  from  the  bottom  boundary  *in  the  first 
column  of  flow  field  cells  along  the  side  face.  This  cell  is  on  the 
downstream  side  of  the  90-degree  expansion  corner  from  the  front  face 
to  the  side  face.  The  BAAL  overpressure  values  are  consistently  less 
than  the  HULL  values,  showing  the  same  behavior  as  that  for  analogous 
flow  field  cells  shown  in  Figures  F.4,  F..5,  and  F..7,  figure  F..12  shows 
the  overpressure  history  of  flow  field  coll  (18,8,4),  the  fourth  cell 
up  from  the  bottom  boundary  in  the  last  column  of  flow  field  cells 
along  the  side  face.  Ibis  cell  is  on  the  upstream  side  of  the  90 -deg roe 
expansion  corner  from  the  side  cnee  to  the  back  face.  As  is  the  case 
for  tho  analogous  flow  field  cells  in  Figures  F..1,  F.2,  and  F..5,  the 
BAAL  overpressure  values  are  consistently  greater  than  the  values  from 
either  of  the  HULL  computations.  Figure  F.13  shows  the  overpressure 
history  of  flow  field  coll  (14,8,3),  the  fifth  cell  along  the  third 
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row  of  flow  field  cells  up  from  the  bottom  boundary  on  the  side  face. 

As  is  the  case  for  the  analogous  flow  field  cell  shown  in  Figure  E.9, 
the  BAAL  and  HULL  calculations  show  good  general  agreement,  as  do  the 
flow  field  cells  away  from  the  edges  of  the  top  and  side  shelter  faces 
(such  as  those  shown  in  Figures  E.14  and  E.10). 

E.  General  View 

Figure  E.15  shows  an  exploded  view  of  the  S-280  shelter,  indicating 
the  pattern  by  which  the  BAAL  computation  of  overpressure  on  the  shelter 
differs  from  that  computed  by  HULL.  The  BAAL  computation  of  overpressure 
is  consistently  higher  than  that  computed  by  HULL  for  flow  field  cells 
on  the  shelter  surfaces  just  upstream  of  a  90-degree  expansion  corner, 
and  consistently  lower  for  flow  field  cells  on  the  shelter  surface  just 
downstream  of  a  90-degree  expansion  corner. 
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Figure  F.2.  Overpressure  for  front  face  cell  (9, 


Figure  E.3.  Overpressure  for  front  face  cell  (9,3,3). 
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Figure  E.4.  Overpressure  for  back  face  cell  (19,3,7). 
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Figure  f.S.  Overpressure  for  back  face  cell  (19,7,3). 


Overpressure  for  back  face  cell  (19,3,3) 
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Figure  F..S.  Overpressure  for  top  face  cell  (18,3,8). 


Figure-  E.9.  Overpressure  for  top  face  cell  (14,7,8). 
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figure. £.11.  Overpressure  for  side  face  eel!  (10,8,4), 
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Figure  E.14.  Overpressure  for  side  face  cell  (14,8,3). 


LIST  OF  SYMBOLS 


c  velocity  of  sound  (m/s) 

h  the  clearing  height  (m) ,  equal  to  the  smaller  of  either  the 
height  of  the  target  or  one-half  of  its  width 

p  static,  or  side-on,  pressure  (kPa) ,  absolute  unless  otherwise 

indicated 

t  time  (s),  where  t  =  0.0  when  the  incident  shock  wa-'e  arrives  at 
the  target  front  face,  as  computed  according  to  theory 

u  particle  velocity  (m/s)  with  respect  to  an  Eulerian  reference 
frame 

D  the  depth  (m)  of  the  target,  measured  in  the  direction  of  travel 
of  the  incident  shock  wave 

1  the  specific  internal  energy  (J/kg) 

M  particle  local  Mach  number 

R  gas  constant  (J/kg-°K) 

T  static,  or  side-on,  temperature  (°K) 

W  wave  velocity  (m/s)  with  respect  to  an  Eulerian  reference  frame 
X  the  direction  of  measure  of  depth  (m) 

Y  the  direction  of  measure  of  width  (m) 

2  the  direction  of  measure  of  height  (m) 

a  a  constant  multiplicative  factor 

B  fp2  -  P]>/2Pi 

Y  the  ratio  of  specific  heats 

n  Courant -Friedrichs- Lcwy  (CFL)  stability  factor 

p  static  density  (kg/m3),  absolute  unless  otherwise  indicated 

Subscripts 

b  the  back  face  of  the  target 
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LIST  OF  SYMBOLS  (Continued) 

clearing  time  when  used  with  the  symbol  "t",  the  "...  time 
required  to  clear  the  front  wall  of  reflection  effects..."2 


drag 

the  drag  phase 

f 

the  front  face  of  the  target 

i 

incident  shock  wave 

max 

a  peak  value 

pmin 

when  used  with  the  symbol  "t", 
pressure  on  a  target  face  to 
having  been  loaded  to  a  peak 

the  time  required  for  the  average 
reach  a  local  minimum  value  after 
value  by  a  shock  wave 

rise 

when  used  with  the  symbol  "t", 
face  to  reach  a  peak  average 
of  the  incident  wave  at  that 

the  time  required  for  a  target 
pressure  after  the  initial  arrival 
face 

top 

the  top  face  of  the  target 

0  stagnation,  or  face-on,  value  after  the  flow  is  isentropically 
brought  to  rest 

1  ambient,  atmospheric,  or  reference  condition,  specifically  that 

region  of  undisturbed  gas  ahead  of  the  oncoming  shock  wave 

2  the  region  behind  the  incident  shock  wave 

3  the  region  representing  the  shock  tube  reservoir 

4  the  region  separated  from  the  shock  tube  roservoir  by  the 

expansion  wave,  and  separated  from  the  shockod  gas  by  the 
contact  discontinuity 

5  the  region  behind  the  reflected  shock  wave 
Superscripts 

*  the  valuo  over  the  reference,  or  ambient,  condition  (e.g.,  p* 
represents  overpressure) 

average  value  over  a  given  face  of  the  target 
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